SOLITON DYNAMICS FOR THE NONLINEAR 
SCHRODINGER EQUATION WITH MAGNETIC FIELD 



MARCO SQUASSINA 



Abstract. The semiclassical regime of a nonlinear focusing Schrodinger equation in pres- 
ence of non-constant electric and magnetic potentials V, A is studied by taking as initial 
datum the ground state solution of an associated autonomous stationary equation. The 
concentration curve of the solutions is a parameterization of the solutions of the second 
order ordinary equation x = —yV{x) — i x B{x), where _B = V x A is the magnetic field 
of a given magnetic potential A. 
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1. Introduction 

The aim of this paper is the study of the asymptotic behaviour of the solutions of the 
semihnear Schrodinger equation with an external magnetic potential A, 



0,(x,O) = 0o(x), xeR"", 



in the semiclassical regime of e going to zero, by choosing a suitable class of initial data 0o 
which is related to the (unique) ground state solution r of an associated elliptic problem. 
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We will show that the evolution (pelt) remains close to r, in a suitable sense (and with 
an explicit convergence rate), locally uniformly in time, in the transition from quantum 
to classical mechanics, namely as e vanishes. This dynamical behaviour is also known as 
soliton dynamics (for a beautiful survey on solitons and their stability, see [54]). Here, i is 
the imaginary unit, e is a small positive parameter playing the role of Planck's constant, 
iV > 1, < p < and V : -> M, A : ^ are an electric and magnetic 

potentials respectively. The magnetic field i? is 5 = V x A in and can be thought (and 
identified) in general dimension 2-form of coefficients (diAj - djAi). The magnetic 
Schrodinger operator which appears in problem (P) formally operates as follows 

(1.1) (tV - A{x)Y(j) = -e^A(j) - —A{x) ■ V0 + |A(x)p0 - Tdiv^.A(x)0, 

and it has been intensively studied in works by J. Avron, I. Herbst and B. Simon around 
1978 (see [4, 5, 6, 45, 49] and references therein). U A = 0, then equation (P) reduces to 

f iedtcj), = -4A0, + V{x)<j)e - X e M^, t > 0, 

\Mx,0) = Mx), a;eM^. 

For equation (1.2), rigorous results about the soliton dynamics were obtained in various 
papers by J.C. Bronski, R.L. Jerrard [9] and S. Keraani [38, 39] via arguments purely based 
on the use of conservation laws satisfied by the equation and by the associated Hamiltonian 
system in built upon the potential V, that is the Newton law 

(1.3) X = -VV{x), x{0) = ^0, x(0) = xo. 

For other achievements about the full dynamics of (1.2), see also [30, 31] (in the framework 
of orbital stability of standing waves) as well as [36, 37] (in the framework of non-integrable 
perturbation of integrable systems). Similar results were investigated in geometric optics 
by a different technique (WKB method), namely writing formally the solution as = 
Ue{x, t)e'^*^^'*^/^, where Ug = Uo + eUi + e^U2 ■ ■ ■ , where 9 and Uj are solutions, respectively, 
of a Hamilton- Jacobi type equation (known as ezA;ona/ equation) and of a system of transport 
equations. 

It is very important to stress that, in the particular case of standing wave solutions 
of (1.2), namely special solutions of (1.2) of the form 

(t)e{x,t) =Ue{x)e--^^\ xeM^, teM+, (^GM), 

where Ue '■ M, there is an enormous literature regarding the semiclassical limit for the 

corresponding elliptic equation 

-^Aue + V{x)ue=\ue\^Pue, xgM^. 

See the recent book [2] by A. Ambrosetti and A. Malchiodi and the references therein. 
Moreover, there are various works admitting the presence of a magnetic potential A, and 
studying the asymptotic profile of the solutions : ^ C to 

i(fV - A{x)Yu, + V{x)Ue = Kl'^PUe, X G M^, 

as e goes to zero (see e.g. [3, 7, 18, 19, 20, 21, 40, 47] and references therein). 

In the special case A = V = 0, the orbital stability for problem (1.2) was proved by T. 
Cazenave and P.L. Lions [16], and by M. Weinstein in [58, 59]. Then, A. Soffer and M. 
Weinstein proved in [50] the asymptotic stability of nonlinear ground states of (1.2). 
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See also the following seminal contributions {in alphabetical order): W.K. Abou Salem 
[1], V. Buslaev and G. Perelman [10, 11], V. Buslaev and C. Sulem [12], J. Frohlich, S. 
Gustafson, L. Jonsson, I.M. Sigal, T.-P. Tsai and H.-T. Yau [25, 26, 27, 28], S. Gustafson 
and M.I. Sigal [32], J. Holmer and Zworski [33, 34], A. Soffer and M. Weinstein [51, 52], 
T.-P. Tsai and H.-T. Yau [55, 56, 57]. See also the references included in these works. 

Now, in presence of a magnetic, some natural questions arise: what is the role played by 
the magnetic field B7 if B plays a significant role, what is the correct Newton equation 
taking the place of (1.3), which characterizes the concentrating curve and drives the dynamic 
in the semiclassical limit? 

As known, a charged particle moving in a magnetic field B feels a sideways force that 
is proportional to the strength of B as well as to its velocity. This force, which is always 
perpendicular to both the velocity of the particle and the magnetic field that created it (a 
particle moving in the direction of B does not experience a force) is known as the Lorentz 
force. Hence, charged particles move in a circle (or more generally, helix) around the field 
lines of B {cyclotron motion). During the motion, B can do no work on a charged particle 
(cannot speed it up or slow it down) although it changes its direction (See figures 1 and 2). 

As a consequence, with the expectation (which arises from the magnetic-free case) that 
in the semiclassical limit the dynamics is governed by the classical Newtonian law, one is 
tempted to say that, in presence of an external magnetic field -B, the right counterpart 
of (1.3) is given by the following Newton equation 

(1.4) i = -W{x) - XX B{x), x{0) = ^0, x{0) = xq, 

agreeing that x has to be interpreted as a matrix operation (H^x) if we are not in M^. 

Only after full completion of the present paper the author discovered that a first result 
(mass and momentum asymptotics) in this direction was obtained, independently, with 
decay assumptions on 5 , by A. Selvitella in [48], showing that, in fact, the above guess is 
the correct interpretation, in the transition process from quantum to classical mechanics. 

On the other hand, in this paper, we improve the result of [48], proving a stronger result, 
which is precisely the one predicted by the WKB method. Roughly speaking, under suitable 
regularity assumptions on V and A, we show that, given the initial position and velocity 
Xo,^o ill I^^) and taking as initial datum for (P) 

(/) Mx) = r(^^)ei[^("«)-("-"o)+"-«ol, X E M^, 

where r E if^(M^) is the unique (up to translation) real ground state solution (bump like) 
of the associated elliptic problem 

{S) - ^Ar + r = |r|2Pr in M^, 

then there exists a family of shift functions 6^ : M"*" — > [0, 27r) such that 

(1.5) Mx,t) = ^(^^Z^je^[^(^-WH--W)+--W+''.W] +^^, X G M^, t > 0, 

locally uniformly in time, as e goes to zero, where we have set ||ci;j||he = 0{e), and being 
\\4>\\b^ = £^"^11 V(/)||^2 + £^~^||</'|li2- In particular, with respect to [48], the convergence 
rate is explicit and of the order e and, as a direct consequence, the concentration center in 
the representation formula (1.5) (expressing the soliton dynamics as guessed by the WKB 
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method) is exactly x{t) (in [48] formula (1.5) is not achievable, being the convergence rate 
undetermined) . 

The magnetic potential A contributes to the phase of the solution, and x{t) is the con- 
centration line (which can be considerably influenced by the presence of B, see the phase 
portraits in figures 1-2). Initial data (/) should also be thought as corresponding to a point 
particle with position xq and velocity ^q. 

In the case where Co = and Xq is a critical point of the potential V, as equation (1.4) 
admits the trivial solution x{t) = Xq and ^(t) = for all t G M"*", formula (1.5) reduces to 



locally uniformly in time, as e goes to zero (see Remark 2.5). In turn, the concentration 
of 0£ is static and takes place around the critical points of V, instead occurring along a 
smooth curve in R^. This is consistent with the literature for the standing wave solutions 
mentioned above. 



The author wish to thank Proff. I.M. Sigal, E.H. Lieb, T.-P. Tsai and S. Keraani for pointing 
out useful information and bibliographic references related to the soliton dynamic of NLS 
equations. The author also thanks Proff. M. Degiovanni and D. Fortunato for feedbacks on 
a preliminary version of the paper. The author was supported by the 2007 MIUR national 
research project entitled: "Variational and Topological Methods in the Study of Nonlinear 
Phenomena" . 



In Section 2, we introduce the functional framework, the tools and the ingredients needed 
to write the statement of the main result of the paper. Theorem 2.4. In Section 3, we 
collect various preparatory results concerning the characterization of the energy levels of 
the problem, in the semiclassical regime. In Section 4, we state the main approximation 
estimates for the solutions. In Section 5, we get two integral identities for the evolution of 
the mass and momentum densities. In Section 6, we obtain the approximation results for 
the mass and momentum densities. In Section 7, we obtain an error estimate. In turn, we 
conclude the proof of the main result of the paper. Theorem 2.4. Finally, In Section 8, we 
summarize the results obtained. 



(1) The imaginary unit is denoted by i. 

(2) The complex conjugate of any number z G C is denoted by z. 

(3) The real part of a number z G C is denoted by ^z. 

(4) The imaginary part of a number z G C is denoted by $5Z. 

(5) For all z,weCit holds ^{zw) = ^{zw). 

(6) For all z,w EC it holds '^{zw) = —'^{zw). 

(7) The symbol M"*" means the positive real line [0, oo). 




Acknowledgements. 



Organization of the paper. 



Main notations. 
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(8) 
(9) 
(10) 



The ordinary inner product between two vectors a, 6 G is denoted by a ■ 6. 
The standard norm of a function u is denoted by ||it||L2- 
The standard L°° norm of a function u is denoted by || 
The symbols dt and dj mean ^ and -rr- respectively. 



(12) 

(13) The symbol C 



02 



a' 



The symbol A means ^ + ■ ■ ■ + • 

i^), for m G N, denotes the space of functions with continuous 
derivatives up to the order m. Sometimes C""(M^) 



\C" 



\a\<m 



is endowed with the norm 
koo < oo. 



(14) 
(15) 



The symbol J f stands for the integral of / over with the Lebesgue measure. 
The symbol C^* denotes the dual space of C^. The norm of a z/ in C^* is 



sup 



(f){x)vdx 



C2 



< 1 



Clearly, C^* contains the space of bounded Radon measures. 

(16) C denotes a generic positive constant, which may vary inside a chain of inequalities. 

(17) The first and second ordinary derivatives of t ^— x{t) are denoted by x and x. 

(18) We use the Landau symbols. In particular 0{e) is a generic function such that the 
limsup of e~^0{e) is finite, as e goes to zero. 



2. Statement of the main result 

2.1. Functional setup and tools. It is quite natural to consider operator (1.1) on the 
Hilbert space H^^e defined by the closure of C^(M^; C) under the scalar product 

where V^u = {Dfu, . . . , Df^u) and Dj = i~^edj — Aj{x), with induced norm 



\u\ 



Ha,, 



tVu ~ A{x)u dx + I V{x)\u\'^dx < oo. 



The dual space of Ha^s is denoted by H'^ ^, while the space if^ ^ is the set of u such that 

'\u\\'ir2 = I|m||?2 + II (-V - A{X)) U\\l2 < OO. 

' " " M (see 



Moreover, to problem (P) it can be naturally associated the functional E : HA,e 
also formula (2.4)) 

Eiu) = ^ J ^\/u-A{x)u^dx + j V{x)\u\^dx- j iMpP+^rfx. 

C) 



Finally, we consider the functional £ : H^(R^; 

1 



1 



2 7 p+1 



associated with (S) 
u\'P+^dx. 



It is a standard fact that the solution r of (S) is the unique (up to translation) solution of 
the following minimization problem 



(2.1) 



S{r) = min{£'(M) : u E H 



klU^}- 
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We also set 

(2.2) m:=\\r\\l2. 

Also, r is radially symmetric and decreasing, belongs to C^(M^) fl if^(R'^), and it decays 
exponentially together with its derivatives up to the order two, that is 

(2.3) |D"r(x)| < Ce-"l"l, xeR^, 
for some a, C > and all < |a| < 2 (see e.g. [8]). 

2.2. Well-posedness and conservation laws. We recall that in [15, Section 9.1] (see also 
[23]), in the particular case = 3 and when the external magnetic field B = (6i, &2) ^3) is 
constant (thus A is linear with respect to x), the (global) well-posedness of problem (P) 
in the (natural) energy space HA,e as well as the H\ ^-regularity of the fiux for H\ ^-initial 
data was investigated (see Proposition 2.2 below) by T. Cazenave, M. Esteban and P.L. 
Lions. Furthermore, in general dimension and for a general (smooth) vector potential 
A, the (local) well-posedness in the energy space Ha^b has been recently studied in [42] by 
L. Michel. We also wish to cite earlier papers by Y. Nakamura and A. Shimomura [43], 
Y. Nakamura [44] as well as the important paper by K. Yajima [60]. In particular, in [43], 
if B has decay assumptions at infinity, the problem is locally solved in the weighted space 
E(2) C H^{R^; C) of functions / in L'^{R^; C) such that ||x°L>^/||i2 < 00 for all a, /5 with 
|q!|, > and < \a + j3\ < 2 (notice that, via the decay (2.3), the initial datum 0o in {I) 
belongs to the space S(2)). Finally, see also [15, Theorems 4.6.5 and 5.5.1] and an abstract 
result. Lemma A.l, in the Appendix of [17], by T. Cazenave and F.B. Weissler. 

In order to prove the main result of the paper, we will assume (among other things) 
that A is globally bounded (together with its higher order derivatives). Clearly with this 
assumption the well-posedness and regularity features for (P) get easier to study. On the 
contrary, if A is unbounded, there are genuine regularity problems and the situation gets 
more involved [22]. 

Definition 2.1. We say that a (sufficiently smooth) vector potential A : — > is 
admissible with respect to problem (P) if the following Proposition 2.2 holds for A. 

Proposition 2.2. [well-posedness statement] Assume that < p < 2/N. Then, for 
every e > and all 0o ^ HA,e, there exists a unique global solution 

<f),eC{R+,HA,e)nC\R+,H'A,e) 

of problem (P) with sup ||0e(t)||_ffA e < Moreover, the mass associated with 4'e{t), 

K{t) = ^ j mt)\'dx, 

as well as the total energy E,{t) = e-^E{(j),{t)) associated with (P) 



dx 



(2.4) E^it) = ^J If V0,(t) - Aix)Mt) 

+ ^/ yi^)\Mt)\'dx-j-^-^ J \Mt)\''^'dx, 

are conserved in time, namely 

Kit) = KiO), Ee{t) = Ee{0), for all t e M+. 
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Remark 2.3. From Proposition 2.2, due to the choice of the initial data (/), the mass Meit) 
is also independent of e. Indeed, via the mass conservation and formula (2.2), 



(2.5) 



1 



A4(t)=A4(0) = ^ / \ct>,{xM dx = — 



1 



X — Xq 

e 



dx 



lL2 



for all e > and t G IR+. 

2.3. The driving Newtonian equation. Given the initial data Xo,^o ^ ■, we consider 

x(t),e(t):M+- 



being the (unique) global (under the regularity assumptions on V and A indicated below) 
solution of the first order differential system 



(2.6) 



i(t) = ew, 

i{t) = -VV{x{t))-m^B{x{t)), 
x(0) = xo, 

namely the second order ODE (1.4). Notice that, setting 



(2.7) 



'H{t) = i^m? + V{x{t)\ te 



7i is a first integral associated with (2.6), namely 

n{t) = H{0), for all t G M+. 

In general dimension A^, this follows by the elementary observation that, as M^{x) is a 
skew-symmetric matrix, we have ^{t) ■ M^{x{t))C,{t) = for all t G M"*". 



2.4. The main result. We consider the following assumptions on the external electric and 
magnetic potentials, V and A. 

(V) V G C3(M^) is positive and ||1/||c3 < cx). 

(A) A G C3(M^;M^) with ||v4||c3 < oo and A is admissible (cf. definition 2.1). 
Consider if^(M^; C) equipped with the scaled norm ||0||h., 

||0||L=^^-^||V0||i.+e-^||0||i.. 
The main result of the paper is the following 

Theorem 2.4. Let r be the ground state solution of problem (S) and let (f>s be the family 
of solutions to problem (P) with initial data (/), for some Xo,^o ^ ■ Let (x(t),^(t)) be 
the global solution to system (2.6). Then there exist 6 > and a locally uniformly bounded 
family of maps 6s '■ — > [0, 27r) such that, if \\A\\c2 < 6, then 

(2.8) Ux,t) = r(^Z^)e^[^(^WH--W)+-?W+^'.W] 
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locally uniformly in time, where uJe G Elg and \\uje\\ms = ^i^), o,s e 0. Furthermore, 
without restrictions on \\A\\c2, we have 

(2.9) \^^^a:,t)\=r{^^^^)+u;i{x,t), 

locally uniformly in time, where Cj^ G Big and \Cj\\^^ < 0{e), as e ^ 0. 

Some comments are now in order. 

Remark 2.5. If Xq is a critical point of V and = 0, then the solution of system (2.6) is 
(x(t),^(t)) = (xo,0) for all t G M+. Then, the conclusion of the previous result reads as 

locally uniformly in time, where G and ||ci;e||e^ = 0{e) as e — * 0. In particular, this is 
consistent with the literature of the standing wave solutions of (P) in presence of a magnetic 
potential A (see e.g. [3, 7, 18, 19, 20, 21, 40] and references included). 

Remark 2.6. In the framework of Theorem 2.4, by the exponential decay of r, it holds 

|0e(x,t)| <Ce^"^^ + |^.(x,t)|. 

For an arbitrarily small 6 > 0, the solution (ps of (P) is expected to decay exponentially 
in the set = {x G : \x — x{t)\ > 6 > 0, for all t G M+j faster and faster as e — 0, 
namely (ps rapidly vanishes far from the concentration curve x{t). 

Remark 2.7. A typical situation in is when the external magnetic field B = {bi, 62, ^3) 
is constant. Without loss of generality, up to a rotation, one can assume that B = (0, 0, b) 
for some 6 G M. Hence, the corresponding vector potential is A{x,y,z) = |(— 0). In 
this case, for harmonic external potentials V, namely 

V{Xi,X2, X3) = ^{ujIxI + UJ2XI + C^gXg), Uj G M, 

system (2.6) reduces to 



(2.10) 



'xi{t) = 






X2{t) = 






xsit) = 






' m = 


-ujfxi{t) 




m = 


-Co'|x2(t) 




M^^ = 


-u;|x3(t). 





It is clear that, setting some fixed values of 00 j and choosing some initial data, enlarging 
the value of the third component b of the magnetic field B (say, from to 60), the original 
periodic orbit turns into a more and more helicoidal pattern. See figures 1 and 2. 

Remark 2.8. By complicating some arguments, assumption (V) could be relaxed. For 
instance V can be written as Vi + V2, being ||Vi||(73 < 00 and V2 sufficiently smooth. The 
idea is to use the cut-off function indicated in (3.5), which is nonzero in the ball of 
containing the region where the orbit x{t) is confined (see [39]). 
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Figure 1. Phase portrait of system (2.10) with ui = 1, uj2 = 1.4, = 1.2, 
6 = (left, no magnetic field) and 6 = 5 (right, weak magnetic field). 




Figure 2. Phase portrait of system (2.10) with ui = 1, uj2 = 1-4, cjs = 1.2, 
6 = 20 (left) and 6 = 60 (right). The effects of the magnetic field increases. 

Remark 2.9. As remarked in [39], the soliton dynamics behaviour breaks down in the 

critical case p = jj. Indeed, in this case, if we choose xq = = 0, V{x) = and A = 0, 

then the modulus of the solution of problem (P) with initial data (poi^) = is given by 
\(l)e{x,t)\ = (cost)-^/V(^) for all x eR^ and tE [0, f ) (see also [13]). 

3. Preliminary facts 

In this section we collect some preliminary result which will allow us to prove the main 
result. Theorem 2.4. 

3.1. Magnetic momentum. The following vector function is useful to pursue our goals. 

Definition 3.1. We define the momentum of the solution (p^, depending upon the vector 
potential A, as a function p^ : x M"*" M^, by setting 

pf{x,t) := j^Q{^,{x,t){eV<Pe{x,t) - iA{x)(f),{x,t)), xER^,tE M+. 

First we state the following 
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Lemma 3.2. Let (f)^ be the solution to problem (P) corresponding to the initial data (/). 
Then there exists a positive constant C such that 



TV0,(-,t)-A(x)0,(-,t) 



L2 



< Ce 



N 



for all t G M'^ and any e > 0. 



Proof. The total energy -E'e(t) is conserved (see Proposition 2.2) and it can be bounded 
independently oi e (see Lemma 3.5). Then, since V is positive, defining C,£{x) := ^^(ex), it 
follows that, for some positive constant C, 



(3.1) 



-vc.(-,t)-A(£x)c.(-,t) ^^-C||c.(-,t)L.,+. 



|2p+2 ^ ^ 
.. I r 57,4-5 ^ "-^ • 

I ■ ■ - ■ ■ ^2 

By combining the diamagnetic inequality (see [23] for a proof) 

'V 



\V\C\\<\{--A{ex)')c 
with the Gagliardo-Nirenberg inequality, setting 



a.e. m 



pN 
2p+2 



G (0, 1), we obtain 



WCei-Mv^r... < ||C.(-,t)lli/||V|C.(-,t)|||^. < ||C.(-,t)|| 



1-1? 
L2 



V 
i 



A{ex))U;t) 



L2 



By the conservation of mass (see Remark 2.3), we deduce that ||Ce('! ^)IIl2 = J^eit) 
independently of e. Hence, for all e > 0, we get 



m, 



\\U;t)\\%t".<C 



-rVCe{;t)-A{ex)U-,t) 



pN 



L2 



Since pN < 2 by assumption, the assertion readily follows from (3.1) and rescaling. 

We have the following summability property for p^{x,t). 
Lemma 3.3. There exists a positive constant C such that 



□ 



sup 



pf{x, t)dx 



< C. 



Proof. Taking into account the inequality of Lemma 3.2 and the mass conservation law, by 
Holder inequality we get 

1 



(x, t)dx 



< / \p^{x,t)\dx< 



< 



tN/2 



Ue{-,t)\W- 



N 



\<Pe{x,t)\ -W(Pe{x,t)- A{x)(Pe{x,t) 



dx 



-N/2 



-V<Pe{-,t)-A{x)(t)e{-,t) 



L2 



for all t G M^. The assertion follows by taking the supremum over positive times. 



□ 



3.2. Energy levels in the semiclassical limit. Let us recall a useful tool (see e.g. [39, 
Lemma 3.3]), which reveals useful in managing various estimates that follow. 



Lemma 3.4. Assume that g : 



is a function of class C"^ 



\g\\c'^ < oo, and that 



r is the ground state solution of (S). Then, as e goes to zero, it holds 



g{ex + y)r'^{x)dx = / g{y)r^{x)dx + 0{e 
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for every y G M fixed. Moreover, 0{e ) is uniform with respect to a family JF c C (R ) 
which is uniformly bounded, that is sup^gjr ||5'||c2 < oo. 

In the next lemma we compute the value of the energy associated with (P)-(J), in the 
semi classical regime. 

Lemma 3.5. Let be the energy associated with the family cf)^ of solutions to problem (P) 
with initial data (J). Then, for every t G M'*', it holds 

E,{t) = £{r) + mn{t) + 0{e^), 

as e goes to zero. 

Proof. Notice that, for all X G M^, we get 



X 



_l_ J. 

e / 

Hence, it follows that 



1 \ e 

^ ~ ^o A^ifA(xo).(x-xo)+x.go1 



e 



A{x). 



1 



tV - A{x)j [ry )e = 



dx 



-N 



Vr 



X — Xq 



dx 



^ / r-{^)\A{x,)+i,\'dx + ^ I r-[^)\A{x)\'dx 



2 

7/v 



r'[^^)A{x)-{A{xo)+io)dx 



= J \S7r{x)fdx+\A{xo)+^o\'^m + j r'^{x)\A{ex + xo)\'^dx 
-2 r^{x)A{ex + xo) ■ {A{xo) + ^o)dx. 



In view of Lemma 3.4, we have 



r'^{x)\A{ex + xo)\^dx = \A{xo)\'^m + C(e^), 
j r\x)A{ex + Xq) ■ {A{xo) + Co)dx = A{xo) ■ {A{xo) + ^o)m + 0{e^). 



Then, 



1 



tV - A{x))(r 



X 



^'^^^[Aixoyix-xo)+x-^o] 



dx 



j \Vr{x)f dx + \A{xo) + W^m + \A{xo)\'^m - 2A{xo) ■ (A{xo) + ^o)m + Oi^' 



\Vr{x)f dx + ml^ol"^ + 0{e^ 
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It turn, by combining the conservation of energy (see Proposition 2.2) and the conservation 
of the function Ti (see definition (2.7)), we get 



EAi) = EM = B^(r(^^)eil-«">-<"^-"'+'^«°A 



= 2^/I(t^--^W)(K 

+ / V{xq + ex)r'^{x)dx 



^ ^\ „i[A(xo)-(x-xo)+x-5o] 



1 



dx 



p+1 



\r{x)\^P+^dx 



\\/r{x)fdx + J V{xo + exy{x)dx - J \r{x)\'^P+^dx + ^m\^o\'^ + 0{e^ 
S{r) + j V{xq + ex)r'^{x)dx + ]^m\^^ + 0{e'^) 



£{r) + mV{xo) + -m\io\^ + 0{e^ 
£{r)+niH{Q) + 0{e^) 



as e goes to zero. 



□ 



Lemma 3.6. Let 0^ be the family of solutions to problem (P) with initial data (J). Let us 
set, for any e > 0, t E IR+ and x E 



(3.2) 



V'.(x, t) = e-i«W-[^-+-'Wle-'^("W)-" (j),{ex + x{t),t) 



where {x(t),^{t)) is the solution of system (2.6). Then 

£{Mt))=EM- I y(x)^i%^rfx+lm|e(t) + A(x(t))p 



\ j \A{x)\ 



- m + M^it)) ■ / pt{x,t)dx - m + M<t)) ■ / M^) 



■dx 



I A{x)-pj{x,t)dx. 



Proof. By a simple change of variable and Remark 2.3, we have 

1 



(3.3) Ueit)\\h = m^^ + ^{t),t)\\i, - ^iiv^. 



{t)\\i2 = Kit) = m, tE 



Hence the mass of ipeit) is conserved during the evolution. Let 



Pe{x,t) 



rAf- 



-53(0,(x,t)V0,(x,t)), xER^.tE 



SOLITON DYNAMICS FOR NLS EQUATIONS WITH MAGNETIC FIELD 13 

be the magnetic-free momentum. A direct computation yields 



1 



1 1 



p+1 



\Mt)\''^'dx - m + Aix{t)) ■ / p,ix,t)dx 



2e^ 
1 



1 



2 



dx 



\Aix)\'\Mt)\'dx + / A{x) ■ 53(0,(t)V0,(t)) 



^m|e(t) + A(x(t))|2-l^^ / {(PM^'^-'dx - m + A{x{t)) ■ I Ps{x,t)dx. 



Then, taking into account the definition of E^{t), we obtain 

£{Mt)) = Ee{t) - I V{x)^-^^^^dx + + A{xm 

- i^{t) + A{x{t)) ■ j p,{x,t)dx 

^ \A{x)\'^\(f)e{x,t)\^dx + I A{x) ■ Pe{x,t)dx 



|2 



2e^ 
Finally, since 

Peix, t) = pfix, t) + t) |^ 

we obtain the desired conclusion. □ 

Next we introduce two important functionals in the dual space of C^. 

Definition 3.7. Let (p^ be the family of solutions to problem (P) with initial data (/) and 
let pf be the corresponding momentum. For any t G M"*", let us define an element n^(-,t) 
in the dual space o/C^(M^;M^) and an element Ill{-,t) in the dual space o/C^(M^;R) by 
setting 

eC\R^;R^) : j Ul{x,t) ■ ip dx = j ^ ■ pf{x,t)dx - m^{x{t)) ■ ^{t), 

\/ip eC\R^;R) : J Ul{x,t)ip dx = ^ </?^^^^rfx - my?(x(t)), 

and all t G M"*". Here x{t),^{t) denote the components of the solution of system (2.6). 

We recall a property of the functional 6y on C^(]R^) (see [39, Lemma 3.1, 3.2]). 
Lemma 3.8. There exist three positive constants Kq, Ki, K2 such that, for all y,z E 

Ki\y - z\< \\5y - 6^\\c-i* < K2\y - z\, 
provided that \\Sy — 6z\\c^* < Kq. 

For a fixed time Tq > (to be chosen later on), let p be a positive constant defined by 
(3.4) p = Ki snp \x{t)\ + Ko 

[0,To] 



N 
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where x{t) is defined in (2.6), the constants Kq and Ki are defined in Lemma 3.8, and let 
X be a C°°(M^) function such that < x ^ 1 and 



(3.5) x{x) = 1 if < p. 

Let us now set, for all t G M'^ and e > 0, 



x{x) = if \x\ > 2p. 



{^{t) + A{x{t)) ■Iil{x,t)dx, 

A{x) ■ Iil{x,t)dx, 

\A{x)\'^Iil{x,t)dx, 

{iit) + A{x{t)) ■ A{x)Iil{x, t)dx, 

V{x)Ul{x,t)dx, 



■ye{t) := mx[t) — I xx\x) rr dx. 



rN 



where x is as in (3.5). 

On the functions uol, we have the following estimate. 



Lemma 3.9. There exists a positive constant C = C{V,A) such that 
(3.6) 

where the function : 



J2Hit)\<cnM, 



(3.7) n,{t) :-- 

(3.8) p^it) 
Moreover 

as e goes to zero. 



I^lix, t)dx 
A{x) ■ Iil{x,t)dx 



is defined as VL^{t) = VL^{t) + pf{t), where 



+ he{t)\, te 



+ sup 




ilvllc3<l 





t e 



0,(0) = Oie' 



Proof. Estimate (3.6) is a simple and direct consequence of the definition of uji{t), Oe(t), of 
the uniform boundedness of ^{t), A{x{t)), namely \^(t)\ + \A{x(t))\ < C and of the fact that 
lll^llc-s < oo and ||A||c3 < oo. Let us now prove that ^£(0) = 0{e'^), as £ — > 0. Recalling 
that the initial data 0o is r{{x — xo)/e)e^^^^^^^°'^'^^~^°^~^^'^°\ in light of Lemma 3.4, for any 
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e C2(M^;M^) such that ||(^||c2 < 1, we infer 

ip{x) ■ lll{x,0)dx = J ip{x) ■ p^{x,0)dx — m(f{xo) ■ ^0 
1 



- / <^(x)-53(0e(x,O)V0e(x,O)) 



^ / If (x) ■ A{x)\(j)eix, 0)1'^ dx -mif{xo) ■ C.0 



^ / ip{x) ■ y4(x)r^ ( — - — j dx- mip{xo) ■ 



(p{xo + ex) ■ {A{xo) + ^o)r^ix)dx 

(p{xo + ex) ■ A{xo + ex)r'^{x)dx — m(p{xo) ■ 
= rrnp{xo) ■ {A{xo) + ^o) - rnip^xo) ■ A{xo) 
-mip{xo)-^o + 0{e') = 0{e'), 

as e goes to zero. In a similar fashion, for any ip G C^(]R^) with Hv^Hc^ ^ 1? "we get 

Lp{x)lll{x,0)dx = — J Lp{x)\(l)e{x,0)\'^dx — mLp{xo) 

(p{xo + ex)r'^{x)dx — rrup{xo) = O^e"^). 

Finally, as xixo) = 1, we have |7e(0)| = \mxo - /(xq + ey)xixo + ey)r'^{y)dy\ < O^e"^), by 
Lemma 3.4. This concludes the proof of the assertion. □ 

At this stage, we are ready to estimate the energy values S{ipe{t))- 

Lemma 3.10. Let ipi, he the function defined in formula (3.2). Then there exists a positive 
constant C such that 

< £{^e{t)) - £{r) < CQe{t) + 0{e^), 

for all t G IR+ and e > 0. 

Proof. By combining the conclusions of Lemma 3.5 and 3.6, we obtain 



S{^P,{t)) - S{r) = mn{t) - / V{x) ^''''Z dx + + 



_ dx + -r-^^'-' 

m + A{x{t))) . / p^{x,t)dx - m + A{x{t)) . I 



e' 



+ \ I \A{x)\''^-^^^^dx + I Aix)-pfix,t)dx + Oie'), 
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for all t G M"*", as e goes to zero. Notice that 

m + M^t)) ■ J pf{x, t)dx = mm\' + mA{x{t)) ■ ^(t) + ^i(t), 

A{x) ■ pf{x, t)dx = mA{x{t)) ■ ^{t) + u;2{t), 



lA(x)\'^M^dx = m\A{xm' + ^3(t), 



|2 

|2 



+ A{x{t)) ■ = m^{t) ■ A{x{t)) + m\A{x{t))\'' + u^{t), 



V(x)^M^^dx = mVixit)) + uj,{t). 
It follows that 

SiiJeit)) - S{r) = ^m\m' + mV{x{t)) - mV{x{t)) - u,{t) + ^m|e(t) + A(a:(t))|2 

- m|^(t)|2 - mA{x{t)) ■ ^{t) - iui{t) - m^{t) ■ A{x{t)) - m\A{x{t))\'^ - uj^it) 

+ lm\Aixm' + ^ + mAixit)) ■ ^ + ^2(t) 

= -cu,{t) + u,{t) + - u.it) - u;,{t) + 0{e'), 

which concludes the proof in light of inequality (3.6) of Lemma 3.9. □ 

4. The approximation result 

Let us first recall a useful and well-established stability property of ground states. 

Proposition 4.1. There exist two positive constants A and C such that, if ^ & H^{M.^\ C) 
is such that ||$||l2 = ||r||i2, where r is the ground state solution of (S), and 

£"($) -^(r) < A, 

then 

(4.1) inf 11$ - e'V(- + y)\\lr<C (£:($) - S{r)) . 

yGR^,i9G[0,27r) 

Proof See [58, 59]. □ 
Next, in view of the previous preparatory work, we can state the reppresentation result. 

Theorem 4.2. Let (ps be the family of solutions to problem (P) with initial data (/) and let 
ipe be the function defined in formula (3.2). Then there exist Sq > , a time T* > , families 
of uniformly bounded functions 9^ : ffi''" [0, 27r), y^ : IR+ and a positive constant C 

such that 

(4.2) <P,{x,t) = em-+Mt)+Ai.im.-.it))f x - ye{t) \ ^^^^^^^ 



where 

M)\\m^<C^/^Ut) + 0{e) 
for all e G (0, Eq) and t G [0, T;). 
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Proof. Since the function {t ^e{t)} defined in formula (3.6) is continuous, for any fixed 
To > and Eq, (Tq > 0, we can define the time (recall here that ^2(0) = 0(£^) as e ^ 0) 



(4.3) 



T; := sup {t G [0,To] : (^^(s) < ctq, for all s G (0,t)} > 0, 



for all e G (0, ^o)- Therefore, by choosing the numbers ctq and Eq sufficiently small, by virtue 
of Lemma 3.10, we conclude that 

< S{ip,{t)) - S{r) < CVt,{t) + 0{e^) < A, for all e G (0, Eq) and t G [0, T*). 

Since ||?/'e(t)||2,2 = ||r||/^2, we are in the right position to exploit the stability property 
of ground states (Proposition 4.1). Hence, there exist two families of uniformly bounded 
functions e,:R+^ [0, 27r) and : M+ ^ such that 

^-i^(^ty[ex+x{t)]^-iA{x{t)yx ^^^^^ ^ _ ^ieAt)r(^x + y,(t)) ||^^ < Cn,{t) + 0{e^), 

for all e G (0,£:o) and any t G [0,T*). In turn, by rescaling and setting 9s{t) := e9s{t) and 
Veit) ■.= x{t)- eveit), we get 



X - ye{t) 



< cne{t) + o{e^ 



for all e G (0,£:o) and t G [0,T^*), namely inequality (4.2), concluding the proof. 



□ 



5. Mass and momentum identities 

In the following lemma we obtain two important identities satisfied by the equation. Only 
after completion of the present paper, that the author discovered the second identity was 
independently obtained in [48]. For the sake of self-containedness we include our proof, 
which uses the first identity and it is shorter. 

Lemma 5.1. Let 0^ he the solution to problem (P) corresponding to the initial data (/). 
Then we have the identity 

(5.1) ^^^(x,t) = -div.p^(x,t), xGM^tGM+. 

Moreover, for all t G M+, we have the identity 



(5.2) j^(x,t)dx = -j pf{x,t) X B{x)dx - J VV{x)^-^^^^^dx, 

where B = S/ x A is the magnetic field associated with the vector potential A. 

Remark 5.2. The momentum identity (5.2), which plays an important role in our asymp- 
totic analysis, can be thought as an extension of the so called Ehrenfest's theorem in presence 
of a magnetic field B. 

Remark 5.3. It follows from the momentum identity (5.2) that for the nonlinear Schrodinger 
equation with no magnetic field (V x A = in M^) and with a constant electric potential 
(VV = in M^) the momentum t J p^{x,t)dx is a constant of motion. 
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Remark 5.4. Concerning the addenda in the right-end side of (5.2), in the semiclassical 
regime, by the upcoming Lemma 6.1, as e — > 0, 



We will show that ye{t) remains close to x{t), for e small (cf. Lemma 6.3). Hence, from the 
right-hand side of (5.2) the Newton equation (2.6) naturally emerges, ruling the dynamics 
of a particle subjected to an electric force F^. = — VV'(a;(t)) and to a magnetic force Fi, = 
—v{t) X B{x{t)), being v = x the velocity. 

Proof. By the exponential decay of r(x), dir{x) and dfjr{x) given by (2.3) and the fact that 
< oo, the initial data (/) belongs to i/^^. Hence, by the regularity (see Proposi- 
tion 2.2), it follows that ^^(t) belongs to H^{R^; C) n H^ ^ for all t > 0. By the standard 
Calderon-Zygmund inequality ||3f^</)e(t)||i2 < C|| A0e(t) for all t (see e.g. [29, Corollary 
9.10]) and since, again, ||A||c;i < oo, for any z, j = 1, . . . , we get 




dx ~ m^{t) X B{y,{t)) + mVV{y,{t)). 



s'Wd^^MmL^ < C\\e'AUt)\\L2 < C\\ (^V - A{x)YMmL^ + C\\A{x) ■ V0.(t)|U. 

+ C|||A(x)|U(t)||L2+C||div,A(x)0,(t)|U2 

<C\\Mt)\\Hl^ + C\\Mt)\\m<^, 



for alH > 0. Hence </)^(t) G H^{R- 



C), for all t > 0. Set, for j = 1, . . . , N, 



{pf)j{x,t) 



1 



Q{(l)e{x,t){edj(l)e{x,t) - iAj{x)(j)e{x,t)). 



e 



■N 



To prove identity (5.1) notice that, on one hand, we have 



N 



div,p^(x,t) = - Y,dAVe)A^^i) 



= -^—^{djMx^t)i£dj4)e{x,t) -iAj{x)^e{x,t)) 



- ^—Q{(l)e{x,t){ed^j(j)e{x,t) -idjAj{x)(l)eix,t) -iAj{x)dj(j)e{x,t)) 



= ^A(x)-3f?(V0,(x,t)0,(x,t)) 

- -^Q{^e{x,t)A4>,{x,t)) + ■^dw,A{x)\^,{x,t)\\ 



SOLITON DYNAMICS FOR NLS EQUATIONS WITH MAGNETIC FIELD 19 

On the other hand, it follows 

^^(x, t) = t) (^V - t) + Vix)M^, t) - t) t)] ) 



^ -!3(0,(x,t)A0,(x,t)) + ^A{x) ■ 3fJ(4(x,t)V0,(x,t)) 



+ ^div,A(x)|0,(x,t)|2. 



Now, concerning second identity, (5.2), for any j = 1, . . . , A^, it holds 

= 5^-^3(9,0,9,0,) + ei"^$J(49,(9,0,)) - ^^.(^)^ 



£1-^3(9,0,9,0,) + £1-^3(9,(0,9,0,)) - £1-^3(9,0,9,0,) - —A 



1 . 



£^^^^^^^ dt 



2£i-^3(9,0,9,-0,) +£1-^3(9,(0,9,0,)) - ^A,(x)^''^^'' 



£^ dt ■ 

The second term integrates to zero. Moreover, taking into account identity (5.1), we get 
TAj{x) '^^'^^ {x,t)dx = I Aj{x)daYxP^{x,t)dx = — f 'VAj{x)-pf{x,t)dx 



' dt 



y^9iA,(x)3(0,(x,t)9^0,(x,t))(ia; 

N 

+ £-^ / ^Ai{x)diAj{x)\(t)e{x,t)\^dx. 



Concerning the first term in the formula for dt{p^)j^ conjugate the equation, multiply it by 
2£-^9,-0, and take the imaginary part. It follows (summation on repeated i indexes) 

2£i-^3(9,0:9,0,) = -£2-^5R(A0:9,0,) +£-^|A(a;)|23?(0;9,0,) 

+ £i-^div,A(x)3(0,9,-0,) + 2£i-^A(x) ■ 3(V0,9,-0,) 
+ 2£-^l^(x)3f?(0;9,0,) - 2£-^|0,|2^'3?(0:9,0,) 

= -£^-^5^(9, (9,019,0,)) +£^-^9, (^^) 

+ £-^|A(x)|23fJ(0:9,0,) + £i-^div,A(x)3(0,9,0,) 
+ 2£i-^A(x) ■ 3(V0,9,0,) +£-^9,- (\/(a;)|0,|2) 

-£-^9,y(x)|0,p-£-^^9,(|0,p^+^). 
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Notice that the following identity can be easily shown (recall that (pelt) G for all t), 



i=l 



Then, recalling that H = {djAi — diAj)ij and that the flux of {/>£ is in H , we infer that 
= -e-^y J2A{x){d,A{x)-diA,{x)Me\'dx 

+ e'~'' I Y,9M^)■'^{^e^^(t)e)dx-e-'' I d,V{x)\(i),\^dx 
i=l 

N 

1=1 

N 

= -e-'' / Y,Mx){djA{x) - diA,{x))\<j)e\^dx 
^ i=l 

N 

J2idjA{x) - diA,{x)) ■ ^{^A<Pe)dx - / djV{x)\(Pe?dx 
1=1 

{M^p^{x,t))jdx - j djV{x)\(l)e\^dx. 



Taking into account the formal identification of the notation —pf[x,t) x B{x) with the 
matrix operation M^p^{x,t), we obtain the assertion. To see this in the three dimensional 
case, recalling that 

we obtain the skew-symmetric matrix 







^2^1-91^2 dsAi 


- d,A; 




-B3 B2 ' 






diA2 


- d2A, 




93A2 


-d2As 




Bs -Bi 






piAs 


- dsAi ^2^3 - 93A2 









-B2 Bi 




Then, setting p\ 


= {p^)i, it follows that 


















-B3 B2 ' 




pI 




yeB2 - 


pIb, 




M.^{x)p^{x,t) = 


Bs -B, 








pIBs- 


pIBi 


= -pf{x,t) X B{x) 






-B2 5i 








pIBi - 


PIB2 





The proof is now concluded. 



□ 



6. Mass and momentum estimates 
First, we have the following control on the mass and momentum. 
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Lemma 6.1. Let Eq > 0, T* > and yeit) be as in Theorem Then there exists a 

positive constant C such that 



-N 



■dx — m6. 



ye{t) 



(C2)* 



pf^^'^\x,t)dx - mat)6y,it) ^^^^^ < Cn,{t) + Oie' 



for every t G [^,T*) and e G (0,£:o)- 



Proof. For any v G if^(]R^;C), we have |V|f|p = \Wv\'^ — ^^'•j'Jp'-" . Then, if ifj^^x^t) is the 
function introduced in formula (3.2), by Lemma 3.10 it follows that 

< Sm) - S{r) + \j '^^^^y^'^ rfx < C^M + 0(5^), 

for every t G [0,Tg) and e G (0,eo)- Moreover, as || I'j/'el ||l2 = ||r||i2, by (2.1) we have 



(6.1) 



-dx < Cn,{t) + 0{e' 



for every t G [0,Tg*) and e G (0,£:o)- Now, by the definition of ip^ (cf. (3.2)), we get 

l^(4v^.)P 

p{(pe{ex + x{t),t)eS/(pe{ex + x{t),t)) - {^{t) + A{x{t))\(t)e{ex + x{t),t)\'^\ 



|0,(£X + x(t),t)|2 

eNpf-i*))^ex + x{t),t) - m\<Pe{ex + x(t), t)p 

|0,(5X + x(t),t)|2 
2A'' 1^*^ (^X -|- x(t) , t) I |^/.\|2 



|</.,(£X + X(t),t)|2 

-25^e(i)-pf"^*"(£a; + x(t),t). 
Hence, by a change of variable, we reach 



\m\\<Pe{ex + x{t),t)Y 



(6.2) 



-dx = e 



Notice that by simple computations, by combining (6.1) and (6.2), it holds 



(6.3) 



_^/, p^^^^'^\x,t) _ Jpf^^'^\x,t)dx\Ux,t)\ 
|0,(x,t)| m 



m 



J pf''^''^\x,t)dx 



m 



<CQ,{t) + 0{e^), 

for every t G [0,T*) and e G (0,£:o)- To prove the assertion, we estimate p£{t), where 
(6.4) p,(t):= [^(x)^M^^dx-m^P{y,) + f pf<'^\x,t)ij{x) - mamiVe 
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for every function ip of class C such that HV^Hc^ < 1- Taking into account that, by the 
definition of fie(t) (cf. formula (3.7)) and \\A — A{x(t))\\c3 < C, we have 



pf''^'^\x,t)dx - m^{t) < pf{x,t)dx-m^{t) + hpf{x,t)-pf^^'^\x,t))dx 



< Cn,{t) + / {A{x) - A(a;(t)))^^^4^rfx 



-AT 



CQ,{t) + I iA{x) - Aix{t))) [M^^ - m5,^^ ]dx 



we can conclude that 



< 



< 



< 



< 



pM^{t))(^x,t)'4!{x)dx - mi{t)^{ye{t)) 

^W)(x,t)[V^(x) - ^{yemdx\ + mvem / pj'^<'^\x,t)dx - miit) 

^(X)|0,(x,t)|2 



m 



+ 



ip{x 



-N 



pA(x(t))| 



x,t)--^( / p^^''^'^\x,t)dx 



for all e G (0,£:o) and t G [0,T^*). Since / p^^^^^^\x,t)dx is bounded (see Lemma 3.3) and 



dx — milj{yi;{t)) 
l0.(x,t)P 



dx 



m l ' 



dx = 0, 



setting i^{x) := ipi^x) — ipiyeit)), it holds 

Pe{t) < / mx)\ -p^ dx + C / \i/j{x)\ dx 



+ / 1^(^)1 



X, t)dx 



rN 



dx + cn,{t). 



From Young inequality and estimate (6.3), it follows 
(6.5) p,(t) < / [Cmx)\ + \\i,{x)\^]^-^^^^dx 



+ 



2- - - . . ^ 

^N/2 P£ ^ ^ ^\xi t) ]_ 

|0e(x,t)| m 

1 

2 



AT 



X, t)(ix 



I0e(x,t)| 



< / [C|^(a;)| + ^|^(a;)r]^^^dx + Ca(t) + 0(.^) 
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Via inequality < 26^ + 2(a - bf with a = e~'^/'^\(f)^{x,t)\ and b = - ye{t))/e), 



e 



C 



:{x,t) \ - r( 



X - Veit) 



dx 



for all e G (0,£:o) and t G [0,T;*), by Lemma 3.4 (as ip{ye{t)) = 0) and Theorem 4.2. 



□ 



Next, we need to show that the distance between the points Veit) found out in the proof 
of Theorem 4.2 and the trajectory x{t) is controlled by ^^(t), as e goes to zero. 

Remark 6.2. We stress that in the proof of the next Lemma we will choose the value of 
To that was introduced in formula (4.3) inside the definition of T*. 

Lemma 6.3. Let yeit) be as in Theorem 4-2. There exist positive constants Eq, ctq and Tq, 
namely the values introduced in (4.3) in the definition of T* such that, for some positive 
constant C, 

\xit)-y,{t)\<Cns{t) + Oie'), 
for all t G [0, T;) and e G (0, Eq). 

Proof. We first show that there exists a time Tq such that |ye(t)| < p, for every t G [0,Tj*) 
with T* < To, where p is the positive constant introduced in formula (3.4). Let us first 
prove that ||(5j/^{t2) — ^ye(ti)\\c^* < P for all ^i? ^2 G [^^T*). Let ip G C^(]R^) be such that 
lly'llc^ < 1- Hence, taking into account Lemma 3.3 and identity (5.1), we get 



|0.(a:,t2)|' |0.(x,ti)p 



-N 



-N 



ip{x)dx 




1 >9|0e| 
dt 



(x, t)ip{x)dtdx 




—(p{x)diYx pfix, t)dtdx 



'tl 

< I|Vv^||l= 



t2 



'V(p{x) ■ pf{x,t)dxdt 



dt j \pf{x,t)\dx <Cy\\c2\t2-ti\<C\t2-ti\. 



Hence, for all ^1,^2^ [0,T*), it holds 



-N 



dx 



■dx 



< C\t2-ti 



In view of Lemma 6.1, the following inequality holds, 

m\\5y^i^t,) - < CTo + Ca(t) + 0{e^) < C{To + ao) 

Here we choose the value of Tq and then of ctq, Eq so small that 

C(To + do) + 0{e^) < mm{mKo,mKoK^}, 



being Kq and Ki the constants introduced in Lemma 3.8. Hence, \ye{t2) — ye{ti)\ < Kq 
for all ^1,^2 £ [0,T^*), and since yeifS) = Xq, we obtain the desired assertion. We can now 
conclude the proof of this Lemma. The properties of the function x imply 



\x{t)-yM < -|7.(t)| + - 
m m 



J0.(x,t)p 
xx[x) dx - mye{t) 
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In light of the first step of the proof, we have xiUeit)) = 1 for all t G [0, T*) and e G (0, Eq) 
so that exploiting again Lemma 6.1, we conclude that 



\x{t)-y,it)\<Cneit) + C\\xx\\c2 
which yields the assertion. 



dx — m6, 



Vsit) 



C2, 



□ 



Finally, we get a strengthened version of Lemma 6.1. 



Lemma 6.4. Let sq > and T* > be as in Theorem 4-2. Then there exists a positive 
constant C such that 



-N 



dx — m6. 



x{t) 



+ t)dx - me(t)5.(t) 11^2, < + 



for every t G [0,T^*) and e G (0,eo)- particular, if \\A\\c2 is sufficiently small, we have 



(6.6) 



-dx — m5 



x{t) 



C2* 



+ (x, t)dx - me(t)5.(t) 11^2* < + 0{e'), 



for every t G [0, T^*) and e G (0, Sq). 

Proof. Notice that, taking into account Lemma 6.1, Lemma 3.8 and Lemma 6.3, we get 



dx — m6. 



x{t) 



< 



C2, 



rN 



■dx — m6: 



Veit) 



C2* 



for every t G [0,T*) and e G (0,eo)- In turn, we also get 

\\pf{x,t)dx - m^{t)6^^t)\\c2* < \\pfix,t)dx - pf''^'^\x,t)dx\ 



IC2. 



\pf^^'^\x,t)dx - mat)Sy^it)\ 



C2* 



C2, 



< sup 

\\v\\c2<'^ 

= sup 

Il'/'llc2<l 



[Aix) - Aixit))M^: 



[Aix) - Aixit))M^) 



dx 



-N 



■dx — m6r 



x(t) 



dx 



< sup \\{A{x) - A{x{t)))ip{x)\\c2 

\\v\\c2<'^ 

— dx - md^{t) 



\<Pe{xM 



-N 



■dx — m5. 



x{t) 



Q2* 



cn,{t) + 

+ cn,{t) + 0{e^] 



< c 



C2, 



cn,{t) + 0{e^) < cn,{t) + 0{e'^ 



for every t G \^,T*) and e G (0,£o)- This concludes the proof of the first assertion. Taking 
into account the definitions of ^eit) and p^it), inequality (6.6) is just a simple consequence. 

□ 



7. Proof of the main result concluded 
In this section we will conclude the proof of the main result. 
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7.1. The error estimate. We now show that the quantity fie(t), introduced in (3.7), can 
be made small at the order 0{e^), uniformly on finite time intervals, as e ^ 0. 



Lemma 7.1. There exists a positive constant C = C{Tq) such that fie(t) < C(To)£:^, for 
all e G (0,60) CLi^d t G [0,T^*). // in addition we assume that \\A\\c-i < 5 for some 5 > 
small, then there exists a positive constant C = C(To) such that fie(t) < C{To)e'^, for all 
e G (0,eo) andt G [0,T*). 



Proof. Taking into account Lemma 6.4, via identity (5.2) of Lemma 5.1, we obtain 



—Ill(x, t)dx 
dt ' 



dp^ 

—^{x,t)dx - m^{t) 



I 

= j p^{x,t) X B{x)dx + j W{x 
- mW{x{t)) - m^{t) X B{x{t)) 

j pf{x,t) X B{x)dx + j Vr(^)^^^^ 



dx 



dx 



rnW{x)5x(t)dx — m / ^{t) x B{x)5x(t)dx 



< 



pf{x,t) -m^(t)4(t)j X B{x)dx 



+ 

< PI|C3 
+ II^I|C3 



pf{x,t)dx - m^(t)6^(t) 
-dx - mdx{t) 



-N 



C2* 



< Ch,{t) + 0{e'^ 



for all e G (0,eo) and t G [0,T*). Hence, recalling Lemma 3.9, it follows that 



nj(x, t)dx 



< 



Ul{x,0)dx 







n^(x, T)dx 



dr 



(7.1) 



< 0{e^)+C / fi,(r)rfr. 
Jo 
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< 1. Then identity (5.1) and Lemma 6.4 yield 



dt e 

(p{x) div^. t)dx — mV (p{x{t)) ■ ^{t) 
Vip{x) ■pf{x,t)dx - j mVifix) ■ i{t)5^{t)dx 
Vip{x) ■ (^pf{x,t) - m^{t)5x{t)^dx 



dx - mV(p{x{t)) ■ ^{t) 



< \\ip\\cs\\pf{x,t)dx-m^{t)5.,^t)\\^,.^ < Cn,{t) + Oie^), 
for all e G (0,eo) and t G [0,T*). Hence, by Lemma 3.9, it follows that 



(7.2) 



sup 

\\<fi\\c3<'^ 



n^(x, t)ip{x)dx 



< sup 

Il¥'llc3<l 



+ sup 

||(p|lc3<1^0 



Ui{x,0)ip{x)dx 



dr 



< 0{e^) + C \ QeiT)dT, 

Jo 

for all e G (0,£:o) and t G [0,T*). Finally, again via identity (5.1) and Lemma 6.4, 

\%{t)\= m^{t) + j xx{x)div^pf{x,t)dx 

fniif) — I V{xx{x)) ■ pf{x,t)dx 



V{xxix))m^{t)6^(^t) - j V{xxix)) ■ p^ {x,t)dx 

< ||V(xx(x))||c2||p^(x,t)rfx-me(t)5.(i)||^,. <Cfi,(t) + 0(£2 
for all e G (0,£:o) and t G [0,T*). This, recalling Lemma 3.9, yields 



(7.3) 

for all e G (0,£:o) and t G [0,T*). By collecting inequalities (7.1), (7.2) and (7.3), we get 

^e{t) < 0{e^) + C [ n,{T)dT 



for all e G (0,eo) and t G [^,T*). Then, by Gronwall Lemma, we have tle{t) < C(To)e^, 
for all e G (0,£:o) and t G [0,T^*). Finally, recalling the definitions of f2e(t) and p^{t) and 
exploiting again Lemma 6.4 concludes the proof. □ 



We are now ready to conclude the proof of Theorem 2.4. 

Let 5 > be as in Lemma 7.1. Let us prove the first part of Theorem 2.4. 

We recall that the value of To > was fixed in the proof of Lemma 6.3 and it just depends 

on the data of the problem, such as V, A, m, N. Moreover, by virtue of Lemma 7.1 and by 
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the definition of T* (see tlie proof of Theorem 4.2), it follows that T* = Tq for all e G (0, Eq), 
up to further reducing the value of ^o- Hence ^^(t) < C{TQ)e^ for all e G (0, ^o) 
t G [0,To]. Now, by Theorem 4.2 there exist two families of functions 6'^ : IR+ [0, 2n) and 
y^:R+ ^ such that 

^ ' V e / He ' 

for all t G [0,To]. On the other hand, by combining Lemma 6.3 with Lemma 7.1, it follows 
that \x{t) — 2/e(t)| < Ce^, for all t G [0,To] and £ G (0, Sq)- Then, taking into account the 
exponential decay of Vr, we obtain 



X — x{t) 



< C 



for all t G [0,To] and e G (0,5o)- Therefore, Theorem 2.4 holds true on the time interval 
[0,To]. Let us take x{Tq) and ^(Tq) as new initial data in system (2.6) and the function 



[X) 



X - x(To) 



\[A{x{To)y{x-x{To))+x-i{To)\ 



as a new initial data for problem (P). Whence, by the previous step of the proof, the 
approximation result holds on the interval [To,2To], and hence on an arbitrary finite time 
interval [0,T], for T > 0. 

In order to prove the second part of the statement of Theorem 2.4 one can follow the 
argument of [48] (essentially relying on [9]). Based upon the identity 



Av 

i 



+ Vk; 



p^[v) := Q{v(yv -iAv)), 



the energy functional rewrites as 

E,it) = Er{t) + E^{t) + E^{t), 

where we have set 

Ent) -.= ^1 V{x)\Mx,t)\'dx, E^{t) := ^1 



\Psix,tW 



E^{t) :- 



2e^ 



\V\(pe\ix,t)\^dx 



p+le 



N 



\Mx,tW 

\Mx,t)\^P+^dx. 



dx. 



Then, following the steps of the proof of [48, Lemma 3.5] (on the basis of the quantitative 
estimate of the expansion of E^ up to a error of 0{e'^), cf. Lemma 3.5), we get 



(7.4) 



0<E^i\<P,\)-E^ir)<CQ,it) + O{e' 



0<E^it)-'-lIP^(^^'^\ 



m 



In turn, the second inequality easily yields 



(7.5) 



,N/2 



\Mx,t)\ 



m 



tN/2 



dx < C(l,{t) + 0{e^ 
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Once inequalities (7.4)-(7.5) holds true, the assertion can be proved by arguing as before. 
In fact, inequality (7.4) yields 

ms\-r{-^^)\\l<CnM + 0{e% 
for some ys{t) E M^. Instead, inequality (7.5) allows to prove inequality (6.6) of Lemma 6.4. 



8. Conclusions 

We have analyzed the soliton dynamics features of subcritical (with respect to global well- 
posedness) nonlinear Schrodinger equations in the semiclassical regime under the effects of 
an external electromagnetic field, showing that the solutions concentrate along a smooth 
curve x(t) : M"*" — which is a parameterization of a solution of the classical Newton 
equation involving a conservative electric force Fe = —W{x(t)) as well as the contribution 
of the Lorenz force Ff, = —x{t) x B{x{t)), being B = V x A the magnetic field. The main 
results improves the results of [48], a recent contribution that the author discovered after 
completion of the paper. The technique is based upon the use of quantum (mass and energy 
for the PDE (P)) and classical ((2.7) for the ODE (1.4)) conservation laws, on the lines 
of an argument introduced in J. Bronski and R. Jerrard in 2000 in [9] making no use of a 
linearization procedure for the equation. On the other hand, the presence of the magnetic 
field introduces new difficulties that have to be handled. Finally, we wish to stress that our 
results are consistent with the current literature regarding the analysis of particular classes 
solutions, such as the standing waves. 
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